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We perform fits to the available charged and neutral kaon production data in e"*" + — >■ A" + X, 
K = K^andKg, and determine the non-singlet combination of kaon fragmentation functions — 
in a model independent way and without any correlations to the other fragmentation functions. 
Only nuclear isospin invariance is assumed. Working with non-singlets allows us to include the data 
at very low momentum fractions, which have so far been excluded in global fits, and to perform a 
first NNLO fit to fragmentation functions. We find that the kaon non-singlet fragmentation function 
■ at large z is larger than that obtained by the other collaborations from global fit analysis and differs 

^SJ ' significantly at low z. 
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I. INTRODUCTION 



Now that the new generation of high energy experiments with the detection of a final state hadron are taking place, 
further tests of QCD and the Standard Model require an accurate knowledge not only of the parton distribution 
functions (PDFs) and as{Mz), but equally of the fragmentation functions (FFs) D'^ {z , fij) . These quantities describe 
7*^ the transition of a parton i at factorization scale /z/ into a hadron h carrying away a fraction z of the parton's 
Oh momentum or energy in the center-of-mass (cm.) frame. Like as{Mz), PDFs and FFs are important quantities 
^ . because they are universal: according to the factorization theorem, once they are known at some suitably defined 
scale /i/ — /i/Oj they can be calculated at any other scale /i/ and used in any type of process. The most reliable way 
to determine them at a given scale is by fitting to inclusive single hadron production data in which the fraction x of 
■ available momentum or energy in the cm. frame carried away by the hadron is measured. 
^ " While there has recently been quite an extensive study on the PDFs, only recently have the FFs received more 
detailed studies, and it has been recognized that alot of uncertainties appear in their determination. 
The most direct way to determine the FFs is the one-particle inclusive e~^e~ annihilation process: 



e+e-^h + X, h = TT^, K^,p/p... (1) 

m 

[ Here and from now on we use the shorthand to mean either a h'^ or is observed (but not both) in a given 
■ event. However, these processes, being proportional to the square of the effective electroweak coupling of the quark 

' ■ • • ?)+ h~ h'^ 

• • . g, determine only the combinations + = -\- , i.e. they cannot distinguish the quark and anti-quark 
FFs. In addition, in the limit of massless quarks, they cannot distinguish between the down-type quark FFs D'^ and 
, which have the same electroweak couplings. Different assumptions are imposed in order to gain more information 
^ _ about the FFs. In order to achieve separate determination of and D^, the semi-inclusive DIS l + N^l + h + X 
and the one-hadron inclusive production processes pp ^ h + X and pp ^ h + X play an essential role. However, in 
these processes the nucleon structure is involved, which introduces further uncertainties. 

At present several sets of FFs are available in the literature such as Kretzer 0, Kniehl-Kramer-Potter 
(KKP) [3], Hirai-Kumano-Nagai-Sudoh (HKNS) 3, de Florian-Sassot-Stratmann (DSS) @|, Albino-Kniehl-Kramer 
(AKK,AKK08) [1,01, etc. Two points should be noted about them: 1) in the DSS and HKNS analyses, different 
relations, based on theoretical prejudice, between different initial FFs have been imposed, and 2) there is significant 
disagreement between the various parametrizations for some of the FFs. It is not clear how much of this disagreement 
can be attributed to the choices of experimental data used by these collaborations and how much to the choice of 
the assumptions imposed on the initial FFs. In this paper we shall consider the possibility of obtaining information 
about the FFs directly from experiment, without any assumptions. 

Recently, in [8], we suggested a model independent approach to FFs. We showed that using only C-invariance of 
strong interactions, the difference cross sections between particle and anti-particle production are expressed solely in 
terms of non-singlet (NS) combinations of the FFs to any order in perturbative QCD. 
There are a number of benefits when performing fits of NS quantities: 

i) There are no statistical correlations with gluon FFs, which introduce the largest uncertainties. 

ii) In their fij evolution they do not mix with other FFs, so the difference cross sections are independent of the 
other FFs at all scales. 
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iii) The NS components do not contain unresummed soft ghion logarithms (SGLs) at small z-values. This allows 
the use of measurements at much lower values of x than in global fit analyses [3-0| ) which would (hopefully) better 
constrain the NS. This would provide stronger tests relative to global fits on the validity of the leading twist calculations 
at small x, where the effects of higher twists, as well as of quark and hadron masses, should be most pronounced. 

iv) A next-next-to leading order (NNLO) fit of the non-singlet components is possible, because the perturbative 
components in the NS sector, namely the splitting and coefficient functions, are known to NNLO. This is in contrast 
to global fit analyses where only next-to leading order (NLO) calculations of cross sections are possible at present. 

Note that here and further on we use the notation x for the measured fraction of the energy E of the process carried 
away by the observed hadron h, while we use z for the fraction of the (unobservable) energy Ep of the fragmenting 
parton carried by the observed hadron: 



X 



g2 



E''/E, z = E'^/Ep, (2) 



i.e. X is the measured quantity, z is the theoretically QCD-defined quantity, E'^ is the cm. energy of the observed 
hadron. In leading order (LO), neglecting transverse momenta and hadron mass corrections, x and z coincide. 

In ig] a model independent approach for determining NS combinations of FFs was developed. It was shown that if 
both charged and neutral kaons are measured in l + N l + K + X, in pp — >■ K+X or in e+e" K+X, K = K^,K°, 
SU(2) isospin invariance of strong interactions implies that the cross section differences a"^ between the charged and 
neutral kaons: 

da^ = - 2a^° (3) 

always determines, without any assumptions about PDFs and FFs, the non-sing let ■ 

In this paper we apply the model independent approach of [8j to the available data on and production in 

e"'"e~ annihilation and determine the kaon non singlet — . This allows us for the first time i) to extract 

-D^ without any assumptions about the unfavoured FFs, commonly used in global fit analysis, ii) to extract 
-Df without any correlations to other FFs, and especially to , iii) to determine — in a larger 
region than in global fits by using all available data, that is typically in the region x > 0.001, and iii) to perform a 
first NNLO extraction of the FFs. Including the small x data should also improve the precision of the FFs at large 
z since, via the convolution in eq. ^ below, all z values in the range x < z < 1 contribute, iv) to perform a first 
phenomenological test of recent NNLO calculations and v) to test, at lower x values than before, the incorporation 
of hadron mass according the procedure of Ref. which becomes more important as x decreases. 

The rest of the paper is organized as follows. In section |ll] we describe our approach to charged and neutral kaon 
production. We show how SU(2) invariance allows to single out the NS combination of the kaon FFs, and our basic 
formula for e"'"e~-kaon production is presented. In section IIIII we describe our method of analysis and justify the 
choice of the parametrizations used. The results of our fits and the comparison with those obtained from global 
fits are discussed in section IIVI The results are summarized in section |Vl Appendix |A] outlines our approach for 
calculating the Mellin transform of harmonic poly logarithms, which is necessary for the NNLO calculations. 



II. OUR FORMALISM 



In this section we describe our approach for extracting the kaon non singlet and contrast it to that in global fits. 
In general, the factorization theorem implies that any inclusive hadron production cross section can be written as 

da^x^ED^Y.jy^^"' (?'^-^/) + ^ 

where Eg is the energy scale of the process, dcr* is the process dependent partonic level cross section for the inclusive 
production of a parton i, determined fully in terms of perturbatively calculable coefficient functions, electroweak 
factors, and of the PDFs for any initial state hadrons, /i/ is the factorization scale, and p > 1. Note that though 
formally dcr* is independent of the renormalization scale ^ that appears as the argument of the running coupling 
flg = as/(27r), it depends on it when calculated in perturbation theory, further we assume /i^ ~ /i^ as usually done. 

In LO the measurable quantity x and the QCD variable z usually coincide because dcr* ( x/z. El, /i^ j cx 5{z ~ x). 
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Although the z dependence of the fragmentation functions D'^Zjiij) is not calculated perturbatively, QCD deter- 
mines perturbatively, via the DGLAP evolution equations, their /i^-dependence: 

-^D':iz,,}) =j:f ^P, D'liz\,}) (5) 



din 



where Pij{z,as) are the perturbatively calculable splitting functions. In addition, the DGLAP equations allow a 
choice of = 0{Es) which prevents the large logarithms \og{Es/ fif) from spoiling the accuracy of the perturbative 
calculations of da^. Flavour and charge conjugation symmetry of QCD allow to combine the quark EEs and quark 
coefficient functions into singlets and non singlets, whose advantage is that they do not mix in their evolution. In this 
paper we shall deal with non-singlets. 

In any kaon-production process, if in addition to the charged ii'^-kaons also the neutral ii'^-kaons are measured, 
no new EEs above those used for are introduced in the cross section. This is a consequence of SU(2) invariance 
of the strong interactions, which relates neutral and charged kaon EEs: 

- -D^^ (6) 

Then for the difference cross section da'^ , eq. ([3]), we obtain the simple expression: 



da'^ix,E!) = J'dzida- ~ da'') {-^,Elf,}) {D^^ - D^K^,)^}) 



(7) 



i.e. in any inclusive hadron production process da"^ always depends only on one NS combination of EEs, namely 
Du - ■ This result relies only on SU(2) invariance for the kaons, eq. dl]), and does not involve any other 
assumptions about PDEs or EEs. It holds in any order in QCD. The explicit expressions for da'^ in e^e~ , SIDIS 
and pp scattering were given in Ref . Q . 

In this paper, we focus on the most precisely measured and most accurately calculated processes 

e+e- ^ {-^,Z) ^ K + X, K = K^,KI, (8) 

for which eq.([7]) reads: 

= A^,ao(s) j\z{el- el)(s)C, -1, a.l/i^)^ {D^^ - Df){z, n)), (9) 

where -/s is the cm. energy of the process, x = 2Eh/y/s, ctq = 47rQ;g,„/s is the Born level cross section for the process 
e+e" — > fJ^^n^ , Nc is the number of colours, and e^(s) is the square of the effective electroweak charge of the quark q: 

eq\s)^el^2eqVeVg^ehz + {vl + al) [{vqf + {aqf] \hz\\ (10) 

with hz = [s/{s — + irnzTz)]/ sin^ 26y\r, iq the charge of the quark q in units of the proton charge, and 

Ve = -l/2 + 2sin^6'vi/, = -1/2, 
q ^ II- 2e, sin2 9w, a, = /|, 1^ = 1/2, /| = -1/2. (11) 



V, 
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We set /i^ = /cs, k — 1, 1/4 and 4 to estimate the theoretical error, i.e. we consider three different choices for ///: 
/i/ — \fsl2\ ^fs and 2-y/s. The energy fraction z is given by z = 2{P^ .q) / q^ = Eh/Ep, Cq is the flavour independent 
perturbatively calculated quark coefficient function: 

Cqiz, fij/s, asiii))) = 5{l - z) + as{ii})&p{z, tx)/s) + 0{al). (12) 
Eq. © is our basic formula which we shall use in our fit to determine (D^* — D^^). 

In our analysis we shall use all available and K% production data presented by the different collaborations 
TASSO [TO]— Ip, MARK II [l^], TPC [H, HRS [l|, CELLO jll, TOPAZ tl6|], ALEPH J7.], DELPHI 
OPAL [il— [iai and SLD [2l| at different values of s. 

Experimental data for hadron production ([S]) are commonly presented as normalized to the total hadron cross 
section utot — cq 'Yliq ^q- From eq.® it is clear that the sensitivity of cr^^, {s)/<7tot to {D!^ — ) is determined 
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by the s-dependence of (e^ — e^)(s)/ e^. In Fig. [T]the quantities e^/ J^q ^^'^ ^d/ ^q shown as functions 
of ^/s, which demonstrates that the biggest contribution would come from data away from the intersections with 
the y/s-axis and the region between them, namely away from 80 < y/s < 110 GeV, i.e. most important for our 
studies would be data for which y/s < 60 GeV. It is unfortunate that at the Z-pole y/s ~ 91, 2 GeV, where the most 
precise and abundant data exist, the kaon cross section difference normalized to atot is an extremely small quantity: 



-0.081. 
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FIG. 1: The normalized electroweak charges ej/(e^ + 2e^) (full line) and e^/(e^ + 2e^) (dashed line) as functions of y^. 



For a large part of the e^e^ reaction data for kaon production, the primary quark (i.e. the quark at the electroweak 
vertex) is "tagged" . Experimentally, various techniques are used to achieve this and we refer the reader to the 
various experimental papers (but see in particular Refs. [isj and j20j). For our calculations, we simply neglect the 
contributions from all processes except those for which the primary quark is tagged. Since this can be achieved by 
setting the electroweak charges of all quarks to zero except the tagged quark, the resulting cross section is scheme 
and scale independent as a physical quantity should be. 

We calculate as(/i^) — /(i)/(/3oi), where L = In^^/Agj-^jj and for / in LO, NLO and NNLO we have: 



/lo = 1, /NLo(i) = 1 - 



I3i In L 



(13) 



/nnlo (L) — /nlo (L) 
The constants /?i are given by (2^ 

/3o = 



11 
"6 



2 In^ i - In L 



\n L 



■Ca - -^Tnuf 



(14) 



(15) 



17 5 



(16) 



1415 
216 



11 



79 



108 



(17) 



where Ca = 3, Cf = | and Tr = 1/2. We fix Aqcd = 226 MeV at both NLO and NNLO and for n/ = 5 This is the 



value of Aqcd obtained in the CTEQ6.6M PDF extraction 123|]. 



III. METHOD OF ANALYSIS OF AND K% DATA SIMULTANEOUSLY 



Our formalism would be easy if we had data on and production at identical values of x and ^/s, with the 
cross sections being normalized in the same way. Then the optimum procedure to constrain the kaon non singlet 
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would be to fit it to the difference between these data. However, apart from the u and d quark tagging probabihties 
from OPAL, this does not hold for the data in general. Data on and Kg are at similar cm. energies ^/s, but 
usually at quite different x values. Therefore we proceed in 4 steps; 

1) We combine the measurements on Kg into seven energy intervals -^5 = 12 - 14.8, 21.5 - 22, 29 - 35, 42.6 - 44, 
58, 91.2 and 183 - 186 GeV and parametrize the x dependence of the cross section da^" for JsT^ -production in each 
interval separately as defined below in eq. ([TOl) . 

2) For each interval of ^/s we calculate da"^ perturbatively, using eq. parametrizing the z dependence of the 
kaon non singlet at a suitable starting scale fif — fj,fQ, as described below in eq. (|20p . 

3) Using and (PO]) we calculate the charged kaon cross section a^" in each energy interval of ^/s through the 
expression: 

daf.t- =da,^,-+2daf4„. (18) 

4) We fit the parameters in da^s (as given in eq. (|19l) ) and the parameters in the kaon non singlet (as given in eq. 
(|20p ) simultaneously to measurements of charged and neutral kaon production in e+e^ reactions. 

We believe that the above approach is the optimum one since it involves performing only one fit. 

Since the perturbative calculation of the cross section difference da^ — 2da^s is free of SGLs, it is expected to 
be valid at much lower values of x than the perturbative calculations of da^ and da^^ separately. Therefore in our 
fits we vary the lowest value of x that the data can take. In global fit analyses the usual minimum bound of a; > 0.1 
- 0.05 was used, but in general we will include data at lower values. 

We parametrize the cross section dcr^s as follows: 

{x,s)^( N{s) + ^E^f).\ (1 _ exp[-c(s) In^ x + d{s) In^ x + e{s) In^ x] (19) 

dx \ lay's J 

where TV, A7V, A, B, c, d and e are seven different parameters that are fitted to the data in each range of ^/s separately. 
The AN{s) / In y/s term is motivated by the dependence of the cross section on ^/s predicted by QCD. In the case 
where data of only one ^/s value exists, namely the data at ^/s — 58 and 91.2 GeV, AA^ is fixed to zero. Otherwise, 
note that no QCD input is used for the calculation of the Kg production data. The motivation behind the choice 
of the parametrization is empirical, although the (1 — x)-^ behaviour at large x and the exp[— cln^ x] at small x also 
follow from resummation in perturbative QCD in these respective regions for ^/s ^ Aqcd- 

Due to flavour symmetry, da^ — 2da^^ vanishes whenever the quark at the electroweak vertex is neither a u 
nor d-quark. Thus, we do not need the s, c and 6-quark tagged data from OPAL that should automatically cancel 
and cannot constrain the kaon non singlet in our approach. However, we shall use the light-quark tagged data, that 
contain the u and d-quarks. We can parametrize these data directly, but instead we parametrize c and b quark tagged 
data as in eq. P^ . and calculate the light quark tagged cross section as the difference between the untagged cross 
section and the sum of the c and b quark tagged cross sections. By including all available heavy quark tagged data in 
this way, we hope to improve our calculation of the light quark tagged data. Thus we have 9 parametrized functions 
in X to describe all the K^ data: seven parametrizations for the untagged data in each yi-energy interval and 2 for 
the c and h quark tagged cross sections at -/i = 91.2 GeV. 

For the calculation of do"^ using eq. © , we require a parametrization for the kaon non singlet at a starting scale 
M/ ~ M/o which satisfies the following conditions: It should exhibit the power-like behaviour as z — > 0. Note that 
the resummed double logarithmic contribution to the splitting functions suggests that a Gaussian behaviour in In z 
at small z occurs only for the gluon and singlet FFs [l^l and we do not assume that this behaviour occurs also for the 
non singlet. The FF should also exhibit the behaviour (1 — z)* as z approaches 1. After trying various parametrization 
that were in accordance with the above requirements, we found that the best parametrization, i.e. the one that gave 
a good fit with all parameters well constrained by the data (meaning that the parameters did not become large) , was 

(^f * - Df){z.l^%) = nz'^il - zf + n'z'"'(l - zf. (20) 

This parametrization is effectively the same as the one used in the latest global fits in [HI, 0| , except that a' 7^ a in 
order to allow a larger function space at small z to be available to the non singlet. 

To be clear, our main fit (discussed in subsection IIVB|) which determines the NS D^_^ proceeds as follows. We 
determine D^^^ in a simultaneous fit to K^ and K'^ production data - we fit the K'g production data to eq. (fTOl) 
and we fit the K^ production data to the difference of eq. ([T9| (multiplied by 2) and da"^ : da^ = 2da^= — da'^ , 
where dcr"*' is calculated from — using eq. ([7]). Note that if all the Kg production data were measured 
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at the same x and ^/s values, and defined in the same way, as the production data, there would be no need for 
eq. (fT9|) — we would simply fit the theoretical calculation of da"^ , eq. ([7]), directly to the measurements of da'^ 
at each measured x and ^/s value. We stress that, despite the theoretical discussion immediately following eq. (|19p . 
the motivation for the parameterization in eq. (1191) is mainly empirical - as we will see in subsection IIV Al such a 
parameterization describes all Kg production data well. We note, however, that different parameterizations will exist 
which give an equally good fit to the Kg production data but give slightly different results. Such a "parameterization 
error" should in any case be less than the errors on the parameters due to the errors on the measurements. We also 
note that a single simultaneous fit of all parameters to all data is the statistically correct approach. For example, 
fitting the parameters in eq. ^T9\i to the Kg production data and then, as a separate fit, fitting the parameters in eq. 
(PU)) to the K^ data only would not take into account the fact that the fitted values of the the parameters in eq. 
carry significant experimental errors. 

In our perturbative calculations, we choose ^i/o = GeV, 5 active fiavours u, d, s, c, b, and Aqcd — 226 MeV. We 
perform all calculations in Mellin space since this approach is numerically more efficient than explicitly performing x 
space convolutions such as that in eq. 

The NNLO perturbative components for the cross section difference can be obtained using the results of '25*] for 
the non singlet coefficient functions and the results of [26] for the difference between the spacelike and timelike non 
singlet splitting functions. The former, as well as the spacelike non singlet splitting functions of [27| . are presented in 
Mellin space as a weighted sum of harmonic sums. The latter is presented in x space as a weighted sum of harmonic 
polylogarithms. Our approach for determining the Mellin transform of these harmonic polylogarithms is discussed in 
Appendix [X] 

Because the effect of the observed hadron's mass is expected to be significant at low x, we incorporate the hadron 
mass effects according to the method of Ref. [gj. In this case, the scaling variable x, which in the factorization theorem 
is defined as the ratio of the detected hadron's light cone momentum to the overall process's, must be distinguished 
from the energy and momentum fractions measured in experiment and given by 

XE = 2Eh/y/s and Xp = 2|pt|/\/s (21) 

respectively. We stress that xe and Xp equal x only when hadron mass effects are neglected. Otherwise, they are 
related to x via: 

^-4)> .. = .(l + 4)- (22) 
sx^J \ sx / 

The cross sections da^^ /dx and d(j^= /dx that determine da"^ /dx, eq. (0), which we are calculating and which enters 
the factorization theorem, are related to the measurable ones da/dxp and da/dxE via 9]: 



da ^ ^ 1 da 

dxp ^' 1 + m\/[sx'^{xp)\ dx 

da , , 1 da 



where da stands for either da^^ or da^s . We exploit the fact that different data groups use different definitions for 
"x" in order to obtain the kaon mass, by using the above relations and fitting the mass ruh- We assume the masses 
of the neutral and charged kaons are equal. We note that fitting 

The total number of free parameters in our fits to charged and neutral kaon data is 66, and the total number of 
data points is 730. 

We could choose to fit the data at a subset of ^/s values and then predict the remaining data using the universality 
of the non singlet FF thus obtained. However, we choose to simultaneously fit all the available data in order to 
maximize the constraints on the parameters appearing in eq. (I^H)) . As we will see in section IIVI the simultaneous 
description of all data with the same fitted non singlet FF turns out to be good, in accordance with the universality 
of FFs. 



IV. RESULTS OF THE ANALYSIS 



First we perform a fit only to the available Kg production data in order to ensure that the parametrization in eq. 
(fT9|) is adequate for the Kg data that we will use in our extractions of the kaon non singlet. Then we perform a 
simultaneous fit to both charged and neutral kaon production data in e+e~ reactions. In the latter case we perform 
our analysis to NLO and NNLO in perturbative QCD. 
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A. Analysis of Kg data 

Here we present our results from a fit to Kg data only, using the parametrizations in eq. (jl9p . The average per 
data point, Xdf^ for each data set is presented in Table H] together with details of the data set. Also shown, where 
applicable, is the value at the global minimum of A for each data set, which after multiplication by the normalization 
error is the most likely systematic deviation of the central values (see Ref . Q for a complete discussion) , and which 
should obey |A| < 1 for a reasonable fit. 

TABLE I: Summary of the measurements for inclusive single Kg production in e"'"e~ reactions. The column labeled "Cross 
section" gives the type of cross section measured, up to the normalization and possible non zero width x bins. The column 
labeled data" gives the number of data. The column labeled "Norm. (%)" gives the normalization uncertainty on the 
data as a percentage. The values of A and Xdf from the fit described in the text are also given. In this fit the fitted mass is 
niK = 320 MeV. 



Collaboration 


Cross section 


Tagging 


(GeV) 


# 
data 


Norm. 

(%) 


Xdp 


A 


TASSO flO] 


da s 


untagged 


14.0 


9 


15 


0.4 


-0.6 


TASSO [28] 


da s 


untagged 


14.8 


9 




0.3 


-0.4 


TASSO [28] 


da s 


untagged 


21.5 


6 




0.0 




TASSO [10] 


da^ 


untagged 


22.0 


6 




0.1 


0.1 


HRS [14] 


da s 


untagged 


29 


13 




3.2 




MAKK ii 


aa ^ 


untagged 


29.0 


O 1 




O.e 


O.O 


l-rKj \i-0\ 


da'^s 


untagged 


9Q 


o 
o 




u.o 






da'^s 


1 1 nl" a ffa'pn 


33.3 


9 


15 


0.7 


0.2 


TASSO [10] 


da'<'s 


untagged 


34.0 


15 




1.4 


-0.1 


TASSO [28] 


da'^'s 


untagged 


34.5 


15 




1.3 




CELLO [15] 


da^'s 


untagged 


35 


11 




0.5 




TASSO [28] 


da'^'s 


untagged 


35 


15 




1.3 




TASSO [28] 


da^'^ 


untagged 


42.6 


15 




0.5 




TOPAZ [16] 


da'^'^ 


untagged 


58 


7 




0.1 




ALEPH [30] 


da^'^ 


untagged 


91.2 


30 


2 


0.5 


-2.3 


DELPHI [31] 


da^'^ 


untagged 


91.2 


26 




0.7 




OPAL [19] 


da^'s 


untagged 


91.2 


20 


6 


1.0 


-1.1 


OPAL [20] 


da'^'s 


c tagged 


91.2 


5 




0.6 




OPAL [20] 


da^'s 


b tagged 


91.2 


5 




1.7 




SLD [32] 


da'^'^ 


untagged 


91.2 


17 




1.1 




SLD [32] 


da'^'^ 


1 tagged 


91.2 


17 




0.6 




SLD [32] 


da'^'^ 


c tagged 


91.2 


17 




0.7 




SLD [32] 


da'<'^ 


b tagged 


91.2 


17 




1.5 




DELPHI [33] 


da'^'^ 


untagged 


189 


10 




0.7 




DELPHI [33j 


da^'^ 


untagged 


183 


8 




1.3 












331 




1.1 





In general, as seen from the Table, our parametrization provides a good description of all but the HRS data, where 
the description is poor. At -^/s = 91.2 GeV, the b quark tagged data appears to be slightly inconsistent with the other 
data. The value of |A| for the ALEPH data is high, but the fit to the other data at ^/s = 91.2 GeV in general is good. 
Otherwise, both Xy)F ^'^'^ 1^1 ^ 1 which suggests that the parametrization in eq. ([T9| is sufficient to represent these 
data. 

In the caption of Table U we quote the fitted kaon mass uik = 320 MeV, which is somewhat smaller than the true 
mass of 498 MeV. However, it is not significantly different from the value 343 MeV obtained in global fit analyses in 
Ref. 0, where it was argued that kaon production through complex decay chains may cause a significant difference 
between the true mass and the fitted mass, when only direct parton fragmentation is assumed in the calculations. 



8 



B. Analysis of Tfg and 

Here we present the results from the combined analysis of and data. 

We implement large x resummation in our NLO analysis. We resum both leading and next-to-leading logarithms 
(LL and NLL respectively), which are all the classes of logarithms appearing at this order. As shown in the AKK08 
fit 0, this significantly improves fits to charged kaon data at larg e x. Resummation in the quark cross section (or 
quark coefficient function) is obtained from the method of Ref. [34] and the results for the unfactorized partonic cross 
section in Ref. |35|], while resummation in the evolution is performed according to the method in Ref. [36l] . 

Thus we apply the two most optimum theoretical tools to our calculations, namely the NLO results with resum- 
mation, and the NNLO without resummation. 

The measured inclusive and production cross sections and the obtained x^p values, both in NLO and 
NNLO, are shown in Table HIl In general, with the exception of a few data sets, in particular the h quark tagged cross 
section measurements, the description of the data is rather good. However, the kaon mass, both in NLO and NNLO, 
is significantly lower than the one obtained in the phenomenological description of the data (see Table U) only, i.e. 
without pcrturbative QCD, also it is significantly lower than the value 343 MeV obtained in Ref. Q. 




FIG. 2: The kaon non singlet FF obtained in this paper at NLO with large x resummation (labeled "AC (NLO-|-res)") and 
from the calculations of the same quantity from the HKNS 4], DSS ^ and AKK08 FF sets. 



Our results for - in NLO are shown in Fig. [2j In the same figure the NLO results from global fits of the 
DSS, HKNS and AKK08 sets are presented as well. As seen from the figure, at z > 0.5 there is an agreement in shape 
among the different plots of the NS, but our NS is in general larger in magnitude. However, they differ significantly 
at z < 0.5. The most striking difference is the negative value for the NS at z < 0.4 obtained in our approach, while 

all global fit parametrizations imply a positive — > 0. 

In Tables IhII — |V]we show the values of the parameters for our main fit, NLO + resummation. However, we caution 
the reader that, because we begin our evolution a.t fj,f = \/2 GeV, due to neglect of higher order NNLO terms, the 
uncertainties on these parameters in the NLO calculation may be very large and thus may depend significantly on 
the method used for solving the DGLAP equations. This uncertainty is approximately equal to the size of the NNLO 
terms. 

In order to understand the origin of the negative value of — obtained from the difference cross sections 

a-'^ = a'^ — 2a^° , we make a comparison of the charged and neutral kaon production data at various ^/s in Fig. 
131 Such a direct comparison is possible because, for these data, the cross section measurements happen to be defined 
the same way, i.e. they are differential in the same variable and normalized in the same way, which is not typical 
for the data in general. In general, the description of these data is good. According to the (data-theory) /theory 
plots, the calculation for the Kg production data tends to overshoot the central values of the data, while for 
production the behaviour is the opposite, but this is not significant relative to the experimental errors. For x > 0.3, 
the calculated charged kaon production exceeds the neutral except when y/s = 91.2 GeV. However, below this region 
in X the opposite behaviour is observed, i.e. a'^ < for x < 0.3, for all -y/s except ~ 91.2 GeV, where and 
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FIG. 3: Comparison of charged and neutral kaon production at various cm. energies. 



are very similar. As Fig. [T] shows, the sign of the difference of the effective electroweak couplings of the u and d quark 
flavours, for all s except around the Z-pole (78 < s < 122 GeV^), is positive, i.e. — > for s > 78 GeV^ and 
s < 112 GeV^, i.e. at the cross sections that give the main contribution to the NS in Then, following the simple 
LO approach in which convolutions are replaced by ordinary products, eq.® implies that < at z < 0.3, 

and — Df > at 2: > 0.3. Of course these rough arguments do not take into account experimental errors, 

which are rather big for kaon production, or convolutions etc., however they do help to verify the result qualitatively. 

Our negative result for at low z, though justified by the above arguments on the data on a'-^ , is however 

in contrast to the intuitive interpretation for favoured u-quark and unfavoured d-quark transitions. In addition, our 
result is quite different from the DSS, AKK08 and HKNS results. There could be several reasons for this, as well as 
for the unexpectedly low values for the kaon mass mK(NLO) — 124 MeV and mif (NNLO) = 55 MeV, shown in Table 
ITTl Most probably it is due to the different assumptions in the parametrizations and to inclusion of the small x-data 
in our fit. The DSS and HNKS collaborations use the assumption that all light quark unfavoured FFs are equal: 
D^^ — D^^ — D^^ — D^^ , while no assumptions were used in the AKK08 fit and in the analysis in this paper, 
denoted by AC. The fact that the DSS and HKNS non singlet FF, which can be written as ^ +D^^ - - Df^ , 
is lower than the others for z > 0.4 in Fig. [5] suggests that and Dj^ may be overestimated in this region when the 
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TABLE II: As in Table U but for the fit to both and Kg production data, in which the perturbative components in the 
cross section differences da^ — 2da^s are calculated in NLO and NNLO. The fitted mass of the kaon was mx(NLO) =124 
MeV and mx(NNLO) =55 MeV. 



V^UildUUI cLLiUll 




laggmg 


(GeV) 


# 
data 


Norm. 

(%) 


Xdf 
NLO 




Xdf 
NNLO 


Annlo 


TASSO [37] 


da''" 


untagged 


12 


3 


20 


1.0 


-1.1 


0.9 


-1.0 


TASSO [38] 




untagged 


14 


9 


8.5 


1.0 


-0.1 


0.9 


-0.2 


TASSO [38] 


da''" 


untagged 


22 


10 


6.3 


0.3 


-0.5 


0.3 


-0.6 


HRS [fj 


da''" 


untagged 


29 


7 




1.9 




2.3 




MARKII [12] 


da"" 


untagged 


29 


6 


12 


2.1 


-1.4 


2.6 


-1.4 


TPC [39] 


da"" 


untagged 


29 


29 




1.2 




1.8 




TASSO [37] 


da"" 


untagged 


30 


5 


20 


0.9 


-1.4 


0.9 


-1.4 


TASSO [11] 


da"" 


untagged 


34 


11 


6 


1.5 


-1.0 


1.6 


-0.9 


TASSO [11] 


da"" 


untagged 


44 


4 


6 


0.1 




0.1 




TOPAZ [16] 


da"" 


untagged 


58 


12 




0.7 




0.7 




ALEPH [17, 30] 


da"" 


untagged 


91.2 


29 


3 


1.3 


-0.6 


1.2 


-0.7 


DELPHI [18] 


da"" 


untagged 


91.2 


23 




0.2 




0.2 




DELPHI [18] 


da"" 


1 tagged 


91.2 


23 




0.8 




0.8 




DELPHI [18] 


da"" 


b tagged 


91.2 


23 




0.5 




0.5 




OPAL [40] 


da"" 


untagged 


91.2 


33 




2.3 




2.5 




OPAL [20] 


da"" 


c tagged 


91.2 


5 




4.6 




4.8 




OPAL [20] 


da"" 


b tagged 


91.2 


5 




4.5 




4.4 




SLD [21] 


da"" 


untagged 


91.2 


36 




1.9 




1.5 




SLD [21] 


da"" 


1 tagged 


91.2 


36 




4.7 




4.0 




SLD [21] 


da"" 


c tagged 


91.2 


36 




2.7 




2.4 




SLD [21] 


da"" 


b tagged 


91.2 


36 




4.7 




4.6 




DELPHI [33] 


da"" 


untagged 


189 


8 




5.2 




5.3 




OPAL [20] 


do"" - 2da"s 


u tagged 


91.2 


5 




1.2 




1.4 




OPAL [20] 


da"" - 2da"s 


d tagged 


91.2 


5 




1.0 




1.5 




TASSO [10] 


da"s 


untagged 


14 


9 


15 


0.4 


-0.2 


0.4 


0.0 


TASSO [28] 


da"s 


untagged 


14.8 


9 




0.6 




0.6 




TASSO [28] 


da"s 


untagged 


21.5 


6 




0.1 




0.1 




TASSO [10] 


da"'^ 


untagged 


22 


6 




0.2 


0.2 


0.3 


0.3 


HRS [14] 


da"'^ 


untagged 


29 


13 




2.9 




3.4 




MARK II 


da"'^ 


untagged 


29 


21 


12 


1.2 


1.2 


1.3 


1.4 


TPC [13] 


da"'l 


untagged 


29 


8 




1.8 




2.3 




TASSO [29] 


da"l 


untagged 


33.3 


9 


15 


0.6 


0.3 


0.7 


0.4 


TASSO [10] 


da"'s 


untagged 


34 


15 




1.3 


0.0 


1.3 


0.0 


TASSO [28] 


da"'^ 


untagged 


34.5 


15 




1.3 




1.2 




CELLO [15] 


da"s 


untagged 


35 


11 




0.6 




0.6 




TASSO [28] 


da"s 


untagged 


35 


15 




1.9 




1.9 




TASSO [28] 


da"s 


untagged 


42.6 


15 




0.6 




0.6 




TOPAZ [16] 


da"s 


untagged 


58 


7 




1.1 




1.0 




ALEPH [30] 


da"s 


untagged 


91.2 


30 


2 


1.5 


1.9 


1.4 


1.7 


DELPHI [31] 


da"s 


untagged 


91.2 


26 




3.0 




2.8 




OPAL [19] 


da"s 


untagged 


91.2 


20 


6 


2.3 


0.4 


2.2 


0.3 


OPAL [20] 


da"s 


c tagged 


91.2 


5 




1.2 




1.3 




OPAL [20] 


da"s 


b tagged 


91.2 


5 




12.9 




12.8 




SLD [32] 


da"'^ 


untagged 


91.2 


17 




3.2 




3.0 




SLD [32] 


da"'^ 


1 tagged 


91.2 


17 




0.8 




0.7 




SLD ^ 


da"'^ 


c tagged 


91.2 


17 




1.2 




1.2 




SLD [32] 


da"'^ 


b tagged 


91.2 


17 




5.1 




5.2 




DELPHI [33] 


da"'^ 


untagged 


183 


8 




1.9 




1.9 




DELPHI [33| 


da"'^ 


untagged 


189 


10 




2.7 




2.7 












730 




2.3 




2.2 
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TABLE 111: The fitted values of the parameters for {D^ — ){z, /i/o) parametrized as in eq. (|20|l . from our main fit. 



Parameter 


Value 


n 


-6.25 


a 


-0.11 


b 


3.12 


n' 


11.13 


a' 


0.60 


b' 


3.01 



TABLE IV: The fitted values of the parameters for da^^^_ /dx{x, s) parametrized as in eq. I|19[) in the different energy intervals 
^/s from our main fit: NLO with resummation. 



energy interval in [GeV] 


N 


AN 


A 


B 


c 


d 


c 


12 < < 14.8 


1.58 X 10"^ 


7.90 X 10"^ 


-17.2 


-2.33 


10.1 


-2.99 


-0.357 


21. 5< < 22 


8.43 X 10^ 


-4.79 X 10^ 


14.5 


12.2 


-7.24 


1.45 


8.81 X 10"^ 


29 < ^ < 35 


-0.444 


3.81 


-2.57 


3.49 


1.79 


-0.696 


9.25 X 10"^ 


42.6 < ^/s < 44 


3.55 X 10'' 


-1.32 X lO'' 


4.91 


7.10 


-1.66 


4.08 X 10^ 


-3.21 X 10^2 


^/s = 58 


6.03 


(fixed) 


-1.23 


6.26 


1.11 


-0.415 


-4.93 X 10"^ 


= 91.2 


16.1 


(fixed) 


-4.88 


-0.681 


1.68 


-0.338 


-2.92 X 10"^ 


183 < ^/s < 189 


126. 


-646. 


-1.34 


5.16 


0.492 


-0.156 


-1.69 X 10"^ 



light quark unfavoured FFs are fixed to be equal to one another. Maybe this could explain the similarity of the results 
for * - ZJf * obtained from the DSS and HKNS fits on one hand, and of AKK08 and AC at z > 0.5 on the other 
hand (see Fig. [5]). The AKK08 and HKNS analyses used no data below x < 0.05, the DSS analysis used only data at 
X > 0.1, vifhile we include data as low as a; ~ 0.001. The discrepancy may also be a result of various low x effects not 
accounted for in the calculation, such as dynamical higher twist, quark mass corrections, etc. However, perhaps the 
most likely reason are the large experimental errors on the NS. The FFs of the various collaborations should be the 
same within the error (composed of the theoretical errors and the (unknown) experimental errors propagated from 
the fitted data to the FF). Thus if we assume that the various FFs are consistent, then the spread of FFs in Fig. 
[5] gives some indication of the error on the FF, and shows the error increasing drastically with decreasing z. This 
argument assumes that the (similar) assumptions made on the FFs in the DSS and HKNS fits are correct. In any 
case, these results warrant further investigation into the validity of the standard approach at low x. It is promising, 
however, that it is possible to fit low x data {x ~ 0.001) using fixed order perturbation theory. 

The negative value of the non singlet FF D^* + D^^ — D^^ — D^* at low z contradicts the physical argument 
that favoured FFs are larger than unfavoured FFs. This behaviour alone is not too serious since FFs in general are not 
physical. However, the second Mellin moment of a FF is physical, in the sense of being factorization scheme and 
scale independent, and can be interpreted as the fraction of momentum of the fragmenting parton a that is carried 
away by hadrons of species h. The second moment of the non singlet FF is expected then to be positive, but from 
our fit using calculations to NLO (NNLO) the result is -0.07 (-0.1). Most likely this is a consequence of the large 
experimental errors at low z. However, it could also arise from a breakdown of perturbation theory, or from effects 
not accounted for in the calculation at low x, if such effects turn out significantly large. 

In order to check our negative result for the NS at small z, we performed a fit in which a parameterization of the 
form nz'^{l — z^, instead of that in eq. (PO]) . was used which yielded a positive NS (i.e. n > 0). However, the result 



TABLE V: The fitted values of the parameters for da^^^_ /dx{x,s) parametrized as in eq. (|19|l for = 91.2 GeV from our 
main fit: NLO with resummation for the c and b tagged data 



the data 


N 


AiV 


A 


B 


c 


d 


e 


c tagged 
6-tagged 


2.62 
0.164 


(fixed) 
(fixed) 


-7.52 
-8.65 


2.39 
2.01 


3.97 
2.95 


-1.01 
-0.474 


-9.35 X 10"^ 
-3.08 X 10"^ 
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X^p = 2.4 was obtained, which corresponded to a of about 100 points above that for our main fit. Thus the 
parametriztion in eq. (j20p is much more favoured by the data. We also performed a fit in which the NS was fixed to 
zero, and obtained Xdf = 2.4 again. Thus a positive, as weh as a zero kaon non singlets are both allowed by data as 
a whole, but the fits are much worse. Note that in our analysis we include data at very low x, which are the most 
accurate data and any deviations of the fit from these data immediately results in higher x^- It is the small x data 
that raises with the zero and nz°(l — z)^ parameterizations. 




FIG. 4: The kaon non singlet FF obtained in this paper at NLO with large x-resummation and at NNLO without resummation 
from fits for which k = ijij/s = fJ-^/s =1/4, 1 and 4. For both the left and right plots, the fits for the NS FF are arranged in 
their increasing magnitude ai z = 0.8 as follows: at NNLO with k = 1 and then with k — 1/4, at NLO with fc = 1, at NNLO 
with fc = 4 and at NLO with k — 4 and finally with k — 1/4. Note that the fit at NNLO with k = 1 yields almost a zero NS, 
i.e. the curve coincides with the z-axis. 



In Fig. H we compare the NLO and NNLO fits for D^^ - * at fif = 10 GeV (Fig. H left) and at /x/ = 91.2 
GeV (Fig. m right) for various choices of the factorization scale k = ^"j/s, k = 1, 1/4, 4. As seen from this figure, 
there is an extremely strong dependence on the choice of k and on the chosen perturbative order - NLO or NNLO. 
The quality of the fits for all curves is good, which indicates that the errors of the available data are too big to 
constrain the non-singlet FFs. With very accurate data, the spread of the three NLO FFs and the spread of the three 
NNLO FFs should be less and, assuming that the theoretical error in the calculations is sufficiently larger than the 
experimental error, both should give some indication of the theoretical error. However, since large x resummation has 
been applied only to the NLO calculation and not in NNLO, this analysis cannot test the perturbative convergence 
by comparing NLO and NNLO calculations. Also the low x-data, for the first time included in an analysis, might 
have caused troubles. Fig. 2] implies only upper and lower bounds on the non-sing let: < D!^_^ < 0.2 at z > 0.35. 

Note that in the used method, the uncertainties of the FFs are due almost completely to the experimental errors on 
the data. The theoretical error in the calculation is relatively negligible here. Having the experimental errors as they 
are, one would not get such good fits to the NS FF from cross section (not cross section difference) measurements, even 
with similarly large errors. One of the reasons is that, due to the large x-logarithms in the singlet / gluon evolution 
and in gluon coefficient functions, the small x-data cannot be included in the analysis. All curves are positive for 
z > 0.5, and the theoretical error in this region is arguably less for the NNLO fits. More accurate data is needed to 
better determine the scale variation. 

The non singlet FF from the NNLO fit with fc = 1 is close to zero for z > 0.5, suggesting that the cross section o"^ 
is too. This is a consequence of the fact that the production data and the Kg production data (multiplied by 2) 
are very close (see Fig. [3]). This behaviour is consistent with our finding above, that a good fit can be obtained with 
the non singlet FF fixed to zero (i.e. all parameters in eq. (I^Ul) fixed to zero so that only the parameters in eq. 
are varied in the fit). 

In Fig. [5l we examine the sensitivity of the cross section difference to the mass irih for two different values of y/s, 
taking the mass to be ruh = ttik, rnKf^, and 2mK in eq. ([22]), where mx is the fitted mass. As shown in the figure, 
for = 10 GeV the calculation becomes sensitive to ruh at x < 0.1, while for i/s = 91.2 the sensitivity sets in at 
X < 0.01 where most of the data lie. Because the calculation is very sensitive to hadron mass effects at low ^/s and 



13 



small X, these effects strongly affects our fits. Conversely, precisely because the hadron mass effects are important for 
the data in this region means that they cannot be neglected. However, other low -^5, small x effects, such as higher 
twist and mass effects of resonances from which the kaon has been produced, will also be absorbed into tuk after it 
has been fitted. More accurate data will be needed to determine how important these other effects are. 



-(Cross section difference), Vi=10 GeV, no tagging 



Cross section difference, Vi=91.2 GeV, no tagging 




FIG. 5: The fitted kaon cross section difference at different cm. energies. Also sliown is tlie same quantity but witli the kaon 



mass ruh varied from its fitted result rriK- Note that the left plot shows the negative cross section difference. 
The curves in the left and right plot are negative (and not shown) for x > 0.4 and 0.3 respectively. 



2a' 



In Fig. [6] we show the effect on the fitted NS of increasing the lower bound in x on the data. The result when no 
cut is imposed is similar to the result with a cut of a; > 0.005, implying that data for which x < 0.005 might not 
impose important constraints. The largest change in the fitted NS is from x > 0.005 to x > 0.01. The fact that this 
difference is so large suggests that a new (but valid) minimum in has been found. This implies that the accuracy 
of the data at a; > 0.1 is not enough to form the difference cross sections with the required precision to determine 
— . In particular, note the unphysical divergence of the FF as z 1, which is caused by the negative values 
of the fitted b and b' parameters in eq. (|20p . which further indicates the inability of the data at x > 0.1 to constrain 
the FF at large z. Only including the large amount of precise small x data, which through convolution determines 
the (I?^ — )(z) not only a.t z = x, but also at all z > x, allows to determine the NS in the whole z region. 



V. SUMMARY 

The cross section difference da'^ = da^^ — 2da^°s determines uniquely the NS D^^ — D^^ without any assump- 
tions. We have extracted — from kaon production in e'^e^ K + X, K — K^, Kg and compared our 
results to those from global fit analyses, namely the DSS, the HNKS and the AKK08 parametrizations. In contrast 
to global fits, in our analysis i) data at much lower values of x, as low as ~ 0.001, could be included in the fit because 
of the absence of SGLs in NS perturbative quantities, ii) calculations could also be performed at NNLO and iii) no 
assumptions about unfavoured FFs were imposed. The quality of the fits were high suggesting that perturbative QCD 
is consistent with the data, including the very low x measurements. However, the fitted kaon mass rrix, on which low 
X cross sections depend strongly, was found to be somewhat lower than the value obtained phenomenologically, i.e. 
without perturbative QCD, from the fit to neutral kaon production data only. The obtained values for — 
at small z are negative and considerably different from those obtained from global analyses. Fits performed using 
NNLO calculations gave lower theoretical errors on — , suggesting stability of the perturbation series for 
the most of the available cross section measurements. 

The current measurements of inclusive kaon production are not at the level of accuracy required to obtain a 
competitive extraction of as{Mz) from the considered cross section differences. Also it is not enough to really 
constrain the non-singlet D^_^. As our fits show, the error on the FF is large (see Fig. 4). But our method is a good 
one, which allows for the first NNLO analysis of inclusive hadron production and it will be useful for future studies. 
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NLO+res, M = 10 GeV NLO+res, M=91.2 GeV 




FIG. 6: The kaon non singlet FF obtained in this paper at NLO from fits for which various cuts on the data were imposed. 



Admittedly the experimental errors are large, still one would not get such good fits to cross section (not cross section 
difference) measurements, even with similarly large errors, due to small x-logarithms in singlet / gluon evolution and 
in gluon coefficient functions. 

With more accurate data in the future and in greater number, it would be nice to see if we could continue to 
describe the low x data well. In particular, at present the most accurate data is that for which y/s = 91.2 GeV, which 
is the least sensitive to the kaon non singlet due to the similarity between the u and d quark effective electroweak 
charges at this energy. However, such an extraction may become possible once the accurate measurements of kaon 
production at BaBar [4l[ have been finalized, because these data are at ^/s = 10.54 GeV where the quark electroweak 
charges are very different. These BaBar data could also significantly improve the constraints on D^^ — . 



Appendix A: Harmonic polylogarithms in Mellin space 

In this section we describe our procedure for obtaining the Mellin transform of harmonic sums, defined as 

[x) ^ I dxifmi{xi) I dX2fm2{x2) ■ ■ ■ / dXnfm„{Xn) (Al) 

Jo Jo Jo 

where 

foix)^-, hix) = -^, f_,(^^) = ^ (A2) 

X 1 — X 1 + X 

Their Mellin transforms can be expressed as a weighted sum of harmonic sums, defined for integer values of the Mellin 
space variable n as 

^ (sgn(fci))"^ ^ (sgn(fc2))"- ^ (sgn(fc3))"- 
Sk.MMUn)^}^ 2^ 2^ ... (A3) 

ni = l '^1 n2 = l "2 n3 = l 

which can be continued to complex n according to the procedure in [4^ . 

Such weighted sums can be found recursively '43*] by determining, in Mellin space, the dependence of a harmonic 
polylogarithm on the same one without the leftmost index. By performing the integration for the Mellin transform 
of Hp -j^{x) = Jq dyfp{y)Hjfi{y) by parts, which gives Hp j^{n) = (-ffp,?^(l) - M [xfp{x)Hjfi{x)] {n))/n, and then 
performing the replacements xf±i{x) = =F (1 — f±i{x)), we find that our desired relations are 

~ Ho^^{l)-HM .... 
HqMW = (A4) 
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(n) 



(A5) 



-ff-l,nt(") 



1+a: 



(n) 



(A6) 



where M [H^{x) / {I + x)] (n) and the "+" distribution M [[H^{x) / {I - x)]^] {n) = M [H^{x) / {I ~ x)] in) - 
M[Hif^[x)/{l — x)] (1) can be calculated from Hjji{n) by expanding 1/(1 ± x) as a series in x before performing 
the Mellin transform. The result |4^] is simply that, because X]r=i("Fl)*'5'st(i)/*^ = '5'=pp,st('^) which follows from eq. 
^3p . where p > here and in what follows, each term of the form 



nP 



becomes 

Furthermore, each term of the form 
becomes 



nP 



M 



H-^{x) 



l±x 



(A7) 



in) (A8) 



(-ly 



nP 



(Tl)" [5±p,^(oo) - 5±,,^(n - 1)] = (_l)«M!!) + (^1)" ,^(oo) - 5±,,^(n)] 



M 



H^ix) 



l±x 



(A9) 



(n) 
(AlO) 



Although the Sij^(oo) are singular, they may be treated in a symbolic sense because H:j^{n), M [i/^(2:)/(l + x)] [n) 
and M [[i7^(a;)/(l — x)]^ (n) are all finite when the real part of n is suitably large. To complete this recursive 
procedure, we require the Mellin transforms of the simplest harmonic polylogarithms, which are given by 



H,{n) - H_,{n) = -(-1)"^^ + 'j^^l (-1)") (AH) 



A Mathematica file for implementing this procedure and for producing FORTRAN programs to cal- 
culate numerical values of harmonic polylogarithms anywhere in Mellin space can be obtained from 
http : //www. desy.de/~simon/HarmonicSums. 
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